THE CARMICHAEL NUMBERS UP TO 10 17 



RICHARD G.E. PINCH 



Abstract. We extend our previous computations to show that there are 585355 
Carmichael numbers up to 10 17 . As before, the numbers were generated 
by a back-tracking search for possible prime factorisations together with a 
"large prime variation". We present further statistics on the distribution of 
Carmichael numbers. 

1. Introduction 

A Carmichael number N is a composite number N with the property that for 
every b prime to N we have fe^" 1 = 1 mod N. It follows that a Carmichael number 
TV must be square- free, with at least three prime factors, and that p — 1\N — 1 for 
every prime p dividing N: conversely, any such N must be a Carmichael number. 

For background on Carmichael numbers and details of previous computations we 
refer to our previous paper [3]: in that paper we described the computation of the 
Carmichael numbers up to 10 15 and presented some statistics. These computations 
have since been extended to 10 16 and now to 10 17 , using similar techniques, and 
we present further statistics. 

The complete list of Carmichael numbers up to 10 is available from the author. 

2. Organisation of the search 
We used improved versions of strategies described in 0]. 

The principal search was a depth-first back-tracking search over possible se- 
quences of primes factors pi,...,pa- Put P r = 111=1 P*> Qf = Tli=r+iPi an ^ 
L r = lcm {pi — 1 : i = 1, . . . , r}. We find that Q r must satisfy the congruence 
N = P r Q r = 1 mod L r and so in particular Qd = pd must satisfy a congruence 
modulo Ld-i- further pd — 1 must be a factor of Pj-i — 1- We modified this to 
terminate the search early at some level r if the modulus L r is large enough to limit 
the possible values of Q r , which may then be factorised directly. 

We also employed the variant based on proposition 2 of which determines the 
finitely many possible pairs {pd-i,Pd) from Pd-2- In practice this was useful only 
when d = 3 allowing us to determine the complete list of Carmichael numbers with 
three prime factors up to 10 18 . The count of these is 35586 (not 35585 as incorrectly 
reported in our preliminary work). 

Finally we employed a different search over large values of Pd, in the range 10 7 < 
Pd < 10 9 , using the property that Pd-i = 1 mod (p^— 1) and factorising the numbers 
in this arithmetic progression less than lO 18 /^- 

3. Statistics 

We have shown that there are 585355 Carmichael numbers up to 10 17 , all with 
at most 11 prime factors. We let C(X) denote the number of Carmichael numbers 
less than X and C(d,X) denote the number with exactly d prime factors. Tabled 
gives the values of C(X) and Tabled the values of C(d, X) for X in powers of 10 
up to 10 17 . 
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Table 1. Distribution of Carmichael numbers up to 10 
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Table 2. Values of C(X) and C(d, X) for d < 10 and X in powers 



of 10 up to 10 



We have used the same methods to calculate the smallest Carmichael number 
Sd with d prime factors for d up to 34. The results are given in Table Heuristics 
suggest that a good approximation for logS^ might be 21og2(d — l)log(rf — 1). 
There is some support for this from Table 

In Tableland Figure ^ we tabulate the function k(X), defined by Pomerance, 
Selfridge and Wagstaff .6! by 

C(X) = Xc Xp (-k(X) l0SXlOgl0Sl0SX 



log log X 

They proved that liminf k > 1 and suggested that limsup/c might be 2, although 
they also observed that within the range of their tables k(X) is decreasing: Pomer- 
ance [7],[E] gave a heuristic argument suggesting that Yank = 1. The decrease in k 
is reversed between 10 13 and 10 14 : see Figure^ We find no clear support from our 
computations for any conjecture on a limiting value of k. 
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260849323075371835669784094383812120359260783810157225730623388382401 



Table 3. The smallest Carmichael numbers with d prime factors 



for d up to 34. 



In TableEJwe also give the ratios C(10 n )/C(lO n - 1 ) investigated by Swift [TUj. 
Swift's ratio, again initially decreasing, also increases again before 10 . 

In Tableland Figure[3we see that within the range of our computations C(X) 
is a slowly growing power of X: about X 339 for X around 10 17 . 

In Table we give the number of Carmichael numbers in each class modulo m 
for m = 5, 7, 11 and 12. 

In Tables |H1 and we give the number of Carmichael numbers divisible by primes 
p up to 97. In Table IS1 we count all Carmichael numbers divisible by p: in Table 
we count only those for which p is the smallest prime factor. 

The largest prime factor of a Carmichael number up to 10 17 is p = 223401361, 
dividing n = 99816335969903281 = 17 • 31 • 379 • 2237 • p. We have p = n 0A911 : 
indeed n = p(2p — 1) and n — 1 = (p — l)(2p + 1). 

The largest prime to occur as the smallest prime factor of a Carmichael number 
in this range is 380251, dividing 



90256390764228001 = 380251 • 410761 • 577981. 
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d 




21og2(d-l) log(d-l) 
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2.324869052 
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1.385943139 
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1.296862582 
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1.273113126 
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1.210794211 
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1.187291825 
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1.183842403 
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1.142841324 
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1.115637251 
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1.112255293 
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1.068846695 


18 


1.086833722 


19 


1.067861876 


20 


1.055266652 


21 


1.056211783 


22 


1.041906608 


23 


1.034833298 


24 


1.024202006 


25 


1.026042173 


26 


1.014073506 


27 


1.017122489 


28 


1.010169037 


29 


1.007224804 


30 


1.000945160 


31 


0.984182030 


32 


0.993535534 


33 


0.990337138 


34 


0.984854273 



Table 4. The smallest Carmichael number with d prime factors 
Sd compared to 2 log 2(d — 1) log(<i — 1). 



We note that this number is of the form (25fc + l)(27fc + l)(38fc + 1) with k = 15210. 

We define the index of a Carmichael number N to be the integer i(N) — (N — 
1)/X(N) where A is the Carmichael function, the exponent of the multiplicative 
group. A result of Somer 9 implies that i(N) — > oo as N runs through Carmichael 
numbers. In Table 1101 we list the Carmichael numbers known to have index less 
than 100: the list is complete for numbers up to 10 17 . 

We define the Lehmer index of a number N to be £(N) = (N—l)/cf)(N). Lehmer's 
totient problem asks whether there is a number of integral Lehmer index greater 
than 1 (that is, whether there is a composite N such that 4>(N) divides TV — 1). Such 
a number would have to be a Carmichael number. See Guy |2] §B37 for further 
details. 

The Carmichael numbers up to 10 17 with Lehmer index > 2 are listed in Table 
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Figure 1. k(X) versus X (expressed in powers of 10). 
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k(10 n ) 


C(10™)/C(10 Il - i ) 
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2.93319 
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2.19547 
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2.07632 


2.286 
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1.97946 


2.688 
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1.93388 
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1.90495 
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1.87989 


2.533 


10 


1.86870 


2.396 


11 


1.86421 


2.330 


12 


1.86377 


2.286 


13 


1.86240 


2.339 


14 


1.86293 


2.319 


15 


1.86301 


2.353 


16 


1.86406 


2.335 


17 


1.86472 


2.373 



Table 5. The function k{X) and growth of C{X) for X = 10", 
n < 17. 
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FIGURE 2. C(X) as a power of X (expressed in powers of 10). 



n 


log C (10™) /(n log 10) 
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.21127 
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.24082 
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.27224 


7 


.28874 


8 


.30082 


9 


.31225 


10 


.31895 


11 


.32336 


12 


.32633 


13 


.32962 


14 


.33217 


15 


.33480 


16 


.33700 


17 


.33926 



Table 6. C(X) as a power of X. 
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25.10 y 


10 11 


10 12 


10 i3 


10 i4 


10 ib 


10 lb 


10 1Y 


5 





203 


312 


627 


1330 


2773 


5814 


12200 


25523 
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1652 


2785 


6575 


15755 


37467 


90167 


215713 


520990 
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82 


154 


327 


702 


1484 


3048 


6094 


12912 
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172 


344 


725 


1463 


3059 


6285 


12845 




4 
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182 


368 


767 


1519 


3124 


6391 


13085 


7 





401 


634 


1334 


2774 


5891 


12691 


27550 


60602 
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1885 
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11447 


28001 


69131 


168856 


414999 
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105 


186 


432 


967 


2109 


4599 


10011 


21961 
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152 


232 


496 


1055 


2178 


4707 


10039 


21918 
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129 


211 


450 
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2122 


4592 


9944 


21832 
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138 


222 


454 


1033 


2224 


4777 


10125 


21911 
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142 


235 


462 


1018 


2181 


4715 


10158 


22132 


11 





335 


547 


1324 


3006 


7032 


16563 


38576 


90731 
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640 


1131 


2770 


6786 


16548 


40891 


100071 


248840 
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139 


217 


473 


1068 


2361 


5338 


12054 


27268 
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142 


220 


457 


1045 


2348 


5319 


12186 


27483 
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104 


187 


442 


1026 


2317 


5261 


11917 


27276 
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152 


243 


466 


1066 


2370 


5316 


12194 


27670 
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116 


198 


440 


1061 


2400 


5384 


12155 


27362 




7 


122 


195 


458 


1023 


2223 


5165 


11853 


26960 




8 


129 


222 


475 


1107 


2450 


5449 


12012 


27639 




9 


131 


218 


465 


1042 


2285 


5179 


11835 


27158 




10 


153 


227 


471 


1049 


2372 


5347 


11830 


26968 


12 


1 


2071 


3462 


7969 


18761 


43760 


103428 


243382 


579215 
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2 


2 


5 


5 


7 




5 


20 


32 


64 


124 


228 


448 


805 


1470 




7 


47 


75 


147 


289 


547 


1027 


1906 


3589 




9 


25 


36 


60 


103 


165 


294 


560 


1018 




11 














4 


10 


25 


56 



Table 7. The number of Carmichacl numbers in each class modulo 
m for m — 5, 7, 11 and 12. 
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p 


25.10 y 


10 u 


10 12 


10 ia 


10 i4 


10 ib 


10 ib 


10 1 ' 


3 


25 


36 


61 


105 


167 


299 


565 


1025 


5 


203 


312 


627 


1330 


2773 


5814 


12200 


25523 


7 


401 


634 


1334 


2774 


5891 


12691 


27550 


60602 


11 


335 


547 


1324 


3006 


7032 


16563 


38576 


90731 


13 


483 


807 


1784 


3998 


9045 


20758 


47785 


110320 


17 


293 


489 


1182 


2817 


6640 


16019 


38302 


91783 


19 


372 


608 


1355 


3345 


7797 


18638 


44389 


106273 


23 


113 


207 


507 


1282 


3135 


7716 


18867 


46612 


29 


194 


336 


832 


2094 


5158 


12721 


31110 


76647 


31 


335 


571 


1320 


3086 


7270 


17382 


41440 


99845 


37 


320 


535 


1270 


2926 


6826 


16220 


38647 


92744 


41 


227 


390 


1001 


2418 


5896 


14344 


34759 


84977 


43 


184 


296 


772 


1920 


4663 


11594 


28650 


70904 


47 


53 


80 


199 


492 


1223 


2873 


6810 


17002 


53 


92 


160 


351 


813 


2041 


5143 


12256 


30560 


59 


26 


41 


92 


262 


644 


1611 


3959 


9613 


61 


269 


453 


1075 


2542 


6047 


14429 


34503 


83268 


67 


110 


178 


407 


1063 


2540 


6306 


15295 


37521 


71 


104 


194 


521 


1320 


3351 


8546 


21485 


53857 


73 


198 


348 


849 


2145 


4925 


11929 


29072 


71137 


79 


64 


107 


247 


686 


1728 


4318 


10693 


26766 


83 


14 


24 


56 


137 


340 


838 


1929 


4640 


89 


68 


131 


320 


788 


1951 


4981 


12178 


30737 


97 


123 


193 


495 


1277 


3123 


7594 


18706 


45542 



Table 8. Primes occurring in Carmichacl numbers. 
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p 


25.10 y 


10 n 


lO 1 " 


10 i3 


10 i4 


10 ib 


10 ib 


10 iY 


3 


25 


36 


61 


105 


167 


299 


565 


1025 


5 


202 


309 


624 


1325 


2765 


5797 


12175 


25481 


7 


364 


579 


1218 


2557 


5461 


11874 


25915 


57459 


11 


263 


428 


1071 


2509 


5979 


14397 


33893 


80745 


13 


237 


431 


1058 


2462 


5699 


13514 


32025 


76256 


17 


117 


206 


496 


1318 


3244 


8114 


20206 


50170 


19 


152 


244 


532 


1401 


3358 


8141 


20020 


49413 


23 


37 


78 


207 


535 


1360 


3317 


8195 


20803 


29 


55 


103 


284 


729 


1822 


4659 


11577 


29149 


31 


101 


168 


390 


876 


2116 


5153 


12575 


30667 


37 


60 


95 


219 


551 


1401 


3418 


8594 


21382 


41 


35 


68 


171 


414 


1092 


2736 


6788 


17275 


43 


35 


65 


168 


403 


943 


2308 


5520 


13636 


47 


14 


16 


36 


81 


195 


459 


1135 


2854 


53 


19 


30 


55 


147 


363 


973 


2327 


5842 


59 


2 


4 


11 


43 


100 


272 


618 


1542 


61 


34 


58 


148 


364 


851 


1978 


4722 


11278 


67 


8 


18 


50 


123 


317 


815 


1950 


4843 


71 


15 


25 


66 


161 


389 


979 


2480 


6178 


73 


14 


28 


68 


175 


406 


1015 


2508 


6277 


79 


4 


10 


17 


66 


175 


467 


1163 


2873 


83 


1 


1 


4 


8 


39 


79 


175 


457 


89 


10 


16 


23 


55 


148 


409 


1003 


2523 


97 


10 


20 


50 


106 


261 


606 


1413 


3445 



Table 9. Primes occurring as least prime factor in Carmichacl numbers. 
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i 


N 


factors 


5 


6601 


7-23-41 


7 


561 


3-11-17 


18 


42018333841 


11-47-1049-77477 


18 


55462177 


17-23-83-1709 


18 


8885251441 


11-47-1109-15497 


21 


10585 


5-29-73 


22 


2465 


5-17-29 


23 


1105 


5-13-17 


25 


11921001 


3 • 29 • 263 • 521 


31 


62745 


3-5-47-89 


37 


11972017 


43 • 433 • 643 


37 


67902031 


43 • 271 • 5827 


39 


334153 


19 • 43 • 409 


43 


52633 


7 • 73 • 103 


44 


15841 


7-31-73 


45 


8911 


7-19-67 


47 


2821 


7-13-31 


48 


1729 


7-13-19 


49 


1208361237478669 


53 • 653 • 26479 • 1318579 


50 


4199932801 


29 • 499 • 503 • 577 


52 


206955841 


17- 71 • 277-619 


53 


1271325841 


17 • 31 • 179 • 13477 


54 


4169867689 


13 • 29 • 383 • 28879 


55 


271794601 


13 • 19 • 743 • 1481 


60 


6840001 


7-17-229-251 


61 


1962804565 


5 • 103 • 149 • 25579 


67 


11985924995083901 


29 -101 -1427 -16349 -175403 


67 


410041 


41 • 73 • 137 


70 


162401 


17-41 • 233 


76 


4752717761 


11 • 17-107-173-1373 


81 


35575075809505 


5 • 197 • 223 • 353 • 458807 


82 


1496405933740345 


5 • 47 • 317 • 40253 • 499027 


83 


142159958924185 


5-37- 107-58379- 123017 


83 


158664761899885 


5 • 37 • 107 • 53987 • 148469 


83 


204370370140285 


5 • 37 • 107 • 48677 • 212099 


83 


24831908105124205 


5 • 29 • 719 • 3023 • 78790717 


85 


171189355538562901 


23 • 71 • 983 • 1031 • 103437589 


90 


3778118040573702001 


11 • 47 • 1051 • 67967 • 102302009 


92 


520178982961 


11 • 29 • 131 • 607 • 20507 


94 


12782849065 


5 -7- 269 -317- 4283 


95 


8956911601 


11 -17 -127 -131 -2879 


97 


5472940991761 


199 • 241 • 863 • 132233 


97 


721574219707441 


167 • 241 • 5039 • 3557977 


97 


83565865434172201 


103 • 1993 • 9551 • 42622169 


97 


9729822470631481 


127-409 • 110681 • 1692407 


99 


438253965870337 


43 • 139 • 49409 • 1484009 



Table 10. Carmichacl numbers with small index 



THE CARMICHAEL NUMBERS UP TO 10 17 



11 



I 


N 


factors 


2.00163 


3852971941960065 


3 


5 


23 


89- 


113 


• 1409 • 788129 


2.00388 


655510549443465 


3 


5 


23 


53- 


389 


• 2663 • 34607 


2.00837 


13462627333098945 


3 


5 


23 


53- 


197 


• 8009 • 466649 


2.01001 


26708253318968145 


3 


5 


17 


113 


• 57839 • 16025297 


2.03067 


269040992399565 


3 


5 


23 


29- 


4637 • 5799149 


2.03207 


158353658932305 


3 


5 


17 


89- 


149 


• 563 • 83177 


2.03614 


1817671359979245 


3 


5 


23 


29- 


359 


• 11027-45893 


2.07138 


16057190782234785 


3 


5 


17 


29 • 


269 


• 6089 • 1325663 


2.07294 


75131642415974145 


3 


5 


23 


29- 


53 • 


617-9857-23297 


2.08129 


881715504450705 


3 


5 


17 


47- 


89 • 


113 • 503 • 14543 


2.08381 


31454143858820145 


3 


5 


17 


23- 


2129 • 39293 • 64109 


2.08894 


6128613921672705 


3 


5 


17 


23- 


353 


• 7673 • 385793 


2.09841 


12301576752408945 


3 


5 


23 


29- 


53 • 


113- 197-1042133 


2.11432 


1886616373665 


3 


5 


17 


23- 


83 • 


353 • 10979 


2.11493 


3193231538989185 


3 


5 


17 


23- 


113 


•167 -2927 -9857 


2.13936 


11947816523586945 


3 


5 


17 


23- 


89 • 


113 -233 -617 -1409 



Table 11. Carmichael numbers with Lehmer index greater than 2 
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